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SUMMARY 

The dynamic analysis for earthquake- and vrind-induced response of a long span, cable- 
stayed freeway bridge by NASTRAN in conjunction with post-processors is describud. 
Details of the structural modeling, the input data generation, and numerical results are 
given. The influence of the dynamic analysis on the bridge design Is traced from the 
project initiation to the development of a successful earthquake and wind resistant con- 
figuration. 

INTRODUCTION 


During the summer of 1972, plans were formulated to design and build a nev/ freeway 
bridge in Seattle, Washington, crossing the lower Duwomish waterwoy. This structure, 
called the West Seattle Freeway Bridge, is to provide a four lane highway ond public 
transit connection between the city and the nearby residential and commetcial areo of 
West Seattle. The Duwomish waterway at this locotlon is novigable by large vessels 
and the bridge is required to be both high and long, so as not to interfere with the 
water traffic. For these and esthetic reasons, a cable-stayed design was decided upon. 
Figure 1 illustrates the initial design concept(U. Plonview and elevation view curvo- 
tures are required by the orientations of the connecting freewoy approach strjctures. 

The main foundation, supporting the tower from which the cobles are suspended, is 
located near the edge of the waterway channel, and all foundations ore supported by 
piles driven into the deep, soft, saturated soil at the site. The initial design incor- 
porated a deck structure consisting of a slab supported by girders and a rigid -frame 
type of tower structure, os shown In the figure. 

Because of Seattle's location in an earthquake zone, and becouse bridges such as the 
West Seattle design ore subject to wind -induced oscillations, it was decided by the 
Seattle City Engineering Department to conduct a thorough dyncenic anolysis. The 
Boeing Company was engaged to perform this analysis. 


(I) Configuration and detail design data shown in this paper were provided by the 
firm of Knoerle, Bender, Stone, and Associates, Inc., Consulting Engineers, 
Seattle, Washington; retained by the City of Seattle to perform the engineering 
design for the West Seattle Freeway Bridge project. 
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The initial work plan included analyses of o number of different bridge preliminary de- 
signs, including both concrete and steel deck constructions, in support of the develop- 
ment of a final structural design concept. Finite clement analysis was decided upon. 
The need for parameter studies of the configurations of the deck, pile foundations, and 
the tower and pier structures was forseen. Thus, a large number of analyses, with 
little setup time between them, was onticipated, and a simple computer model was de- 
sired. However, the nature of the deck design, particularly the combination of gir- 
ders v/ith a long span curved in both elevation and planform, suggested that a complex 
structural behavior might occur, requiring a correspondingly careful structural modeling. 
Therefore it was decided to perform two types of analyses, the first quick and simple, 
representing the deck structure as a single curved beam, and the second more detailed, 
representing the structural components of the deck by individual finite elements in the 
computer model. Figures 2 and 3 are computer plots which illustrate these two mo- 
dels. Though Figure 3 is quite crowded with element lines, the individual girder web 
and flange elements can be seen at the right end of the span. A verification of the 
validity of the data computed with the simple model was planned to be obtained by 
a comparison of its modes and frequencies with those of the complex mode. This, in 
addition to arranging the coniputer coding to facilitate convenience in parameter stu- 
dies, constituted the overall work olan for the structural modeling. 

The specific goals of the dynamic analysis were the calculation of earthquake-induced 
stresses in the structure and the calculation of the critical windspeeds at which aero- 
dynamically induced unstable deck oscillations could occur. Predictions of these data 
were made on the sequence of bridge designs which were generated as the project 
developed. It was found that both earthquake- and wind-induced responses are critical 
design conditions, arxi that the initial types of design configurations were not capable 
of withstanding these responses (Reference 1). On the basis of these early evaluations, 
criteria were developed for achieving dynamically sotisfactory designs. Principally, 
these criteria specified the frequencies of the structural vibration modes to avoid 
large earthquake response and specified the deck torsional stiffness and the shape of 
the deck cross-section to avoid wind-induced unstable oscillations. 

The criteria led to a modification of the deck cross-section to the slant-sided, multi- 
cell closed section shown in Figure 4a and to a modification of the t^wer to the wall- 
type configuration shown in Figure 4b. Designs incorporating the features of Figure 4 
ore satisfoctory for both wirnl- and earthquake- induced dynamic response. Currently, 
additional studies are underway to optimize this basic design for earthquake resistance 
by adjusting its vibration mode frequencies to avoid the known frequencies of princi- 
pal earthquake excitation. This work has achieved a significant reduction in the re- 
quired reinforcing steel. The detail design phase of the work will continue to be 
supported by dynamic analysis until the design is finalized, in the foil of 1973. 

This paper discusses those aspects of the work which ore associated with the finite ele- 
ment idealization and the modal anolysis. This work has been done with the NASTRAN 
system, which hos proved to be a highly effective tool in this application. 
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I SELECTION OF NASTRAN 

I 

t The NASTRAN system has several features which are advantageous for thir problem, 
Heading to its choice over other available structural analyzers. Several of these fea- 
tures are mentioned briefly in this section; others are discussed in more detail in later 
^descriptions of the finite element idealizations. The use of combined cylindrical and 
'Cartesian coordinates was helpful in modeling the combined circular arc and straight 
^line geometry of the bridge planform. The ability to specify nonstructurol mass on the 
tCBAR elements was useful since structural data were provided in the form of mass per 
irunning foot on major structural members. The multi-point constraint feature was parti- 
Scularly useful for representing the connectivity between different portions of the bridge 
Istructure. In particular, this was necessary to represent proper connections of the 
»eck structure to the supporting cables, piers, and the tower, and to represent the 
■footing connections to the piles. Multipoint constraints were also used in the more 
icomptex model (called the 3D model) to connect the individual girders to the deck 
■slab. The NASTRAN plotting feature was used to obtain pictorial descriptions of the 
ptructural vibration modes. In the case of seismic analysis this is particularly impor- 
Itant, because seismic re^onse is strongly dependent on both the shape and the direction 
of principal modal motions. Hence, pictorial data permit a quick, ^qualitative assess- 
ment of the likely seismic importance of the structural modes. And Hnally, since the 
seismic and flutter analyses were done by additional processing of the results of the 
modal analysis, a convenient data access system such as the NASTRAN checkpoint/ 
restart tape feature was required. Thus NASTRAN appeared to be particularly well 
suited to the technical requirements of the problem. 

Additional motivation for using NASTRAN was provided by the ovailobility of Boeing's 
input language, SAIL, which has been adapted for NASTRAN input. Most of the bulk 
data were generated automatically by the use of SAIL. Bridge geometric data were 
provided in equation form, which can be coded directly in SAIL's automatic grid point 
generation format. In addition, SAIL has the capability to generate data within 
special (parameter-controlled) subroutines, called external data generators. External 
data generators were used to generate NASTRAN multi -point constraint equations, the 
girder plus slab deck simulation of the 3D model, and the pile foundation simulations. 
These data generation routines are designed such that a set of input porameters con- 
trols the generation of data. By changing a few of these parameters, a complete new 
let of data can be produced, simulating a new design concept. Through such auto- 
matic input generation it was possible to obtain rapid turnaround of onolyses to support 
the design development. 

A separate computer progrom was written to perform the seismic onalysis, using the 
response spectrum analysis method. This analysis requires the mode shopes, frequencies, 
generalized mass, and internal element forces produced in the modal analysis. The 

I NASTRAN checkpoint/restart tope provided access to these data. Since NASTRAN 
lonnally does not checkpoint element forces, a simple Alter wos used to checkpoint 
the element force file OEFl. 
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BEAM (STICK) FINITE ELEMENT MODEL 


Purpose 

To perform dynamic analysis in support of design trade studies, a simple model of the 
bridge structure, easily modified and with reasonably short run time, was set up. 

Called the "stick" model, it uses simple, beam type representations of all structural 
components. Because of the simplicity «. f this model, it was possible to make para* 
metric studies of important porometers, such as the tower stiffness and the earth lateral 
resistance to pile motions, in order to assess the importance of these factors early 
in the program. 

Description of the Model 

There were two basic configurations from which the porometer studies were made: the 
stdel bridge alternate arxi the concrete bridge alterrKite. The models of each of these 
configurotions included the tower, piers (four piers in the steel model, three in the 
concrete model), the footings, piling, and earth springs to represent lateral earth resis- 
tance to pile motions, and the deck itself. The modeling of the piers, footings, piles 
and earth springs was of primary importance in seismic response. The deck modes arKi 
consequently the deck modeling were of primary interest in the flutter analysis. The 
finite element models included the main span portion of the overall bridge structure, 
which is defined by the locations at which the deck bending continuity with the ap- 
proach spans ts fprminated. This arrangement resulted in analyzing the main span plus 
several shorter adjacent *oans, as required for the particular configuration in question. 

Figures 5 and 6 illustrate th.e steel and concrete finite element models. The outward 
appearances of the ^wo models ore alike except that the concrete model has one less 
pier and a slightly modified tower appearance. In reality, however, the deck proper- 
ties obviously change as do the grid point locotions and all mass and stiffness proper- 
ties. 

The deck geometry in the plon view is a straight line for somewhat less than half the 
span and a circular arc for the remaining part. The steel bridge initioMy analyzed is 
1215 feet long artd the concrete bridge is 1040 feet long. The deck describes a poro- 
bolic ore vertically with a peak elevation of 156 feet. The deck structure consists of 
the concrete slab of the roodbed and the integrally constructor' crete or steel sup- 
porting girders. The outermost girders ore of fascia box const' < .* jn. Figure 1C showe 
o typical croM section for the concrete alterr>o*e. The section properties change olortg 
the span at required by the design moments and the applicable design code loodirtg 
conditions. 

The tower supports the deck throi>gh pin supports . The four cable stays attach to the 
top of the tower orto to the deck 175 feet or either side. The piers support the deck 
by roller type supports which permit relotive .ongitudinol motions. The tower and pier 
footings are supported on pile groups which very in size with 906 piles maximum for 
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the tower and 66 piles minimum for one of the piers. 

The piling simulation is shown by Figure 7. Both pile elements and elements represent- 
ing earth lateral stiffness ore employed. Four s'mu'ated piles can provide a correct re- 
presentation of pile group behavior^ and this number was chosen to keep to a minimum 
the number of elements in the model. Each of the four piles ir a simulated pile group 
consists of two ''.BAR elements, and is fixed at the base (called the point of fixity) and 
provided with four springs (oriented parallel and perpendicular to the spon) ro resist mo- 
tion relotive to the surrounding earth. In the figure, for simplicity, earth springs are 
shown on only one pile. The point of fixity is determined from a detailed pile deflec- 
tion anolysis*, and is the uppermost point at which zero pile bending slope occurs to- 
gether with a very small pile deflection. 

The earth spring element properties ore defined by effectively integroting the distributed 
earth lateral stiffriess over a pile length which is considered tributo'y to c particular pile 
grid point. The application of earth lateral resistance at only two points on the pile is 
on approximation of a type customarily mode in discrete element analysis, and would not 
normally be a cause for concern. In the present case, however, because lateral earth 
stiffness was found to be a very Importont parameter, it was desired to verify the ode- 
quacy of the discrete representation. This was accomplished by comparing pile deflec- 
tions computed for the two grid point pile to those obtained for a many grid point - 
many eorth spring representation. The two grid point pile was fouiJ to predict deflec- 
tion within 10% at the top of the pile. This accuracy is suitable for the dynamic on- 
olysis, and further refinement wiihin the framework of linear elastic analysis does not 
appear worthwhile. 

The simulated pile and earth spring stiffnesses in the finite element model are determined 
to provide the actual combined stiffness of the entire pile group. Denoting by 
earth spring stiffness which would be computed for one actual pile within a pile group, 
the following is the spring stiffnesses required in the finite element simulation. 

k, = k, . N . 1 
^model 7 

N is the number of piles in the grcup. The 1/4 foctor distributes the total group earth 
loterol stiffness to the four simulotev^ piles, in addition, the piles ore located withii; 
the footing oreo (Figure 7) such thot the moments of inertia of le simulated pile oreas 
about the footing longitudinal and transverse oxes match those of the actual pile group. 
This provides simulation of pile group bending stiffnesses. The use of N in the k, for- 
mula would appeor to presume that all piles in the group sustain equal lateral locras 
from the earth. Since this is known to be untrue, an ad{ustment was mode in the k^ 
earth property to occount for group pile oction. The combination of group action and 
vibrotory behavior in the eorthquaK> was occounteJ for by taking k^ to be one-sixth the 
stotic, single pile value. This adjustment is based on reported research on group pile 

* Point of ^xlty ond eorth loterol stiffness doto were provided Ky the firm of Shannon 
orxf Wilson, Soil Mechonics and Foundation Engineers, Seattle, Woshington, 




action . However, In application to a particular pile group, such a factor is neces- 
sarily arbitrary, and It was felt necessary to evaluate the sensitivity of the structural 
behavior to variations in kj^. To accomplish this, a number of computer runs were made 
with widely varying earth spring stiffnesses. It was found that modal and earthquake 
response data are very sensitive to changes in earth spring stiffness for the case of rela- 
tively soft springs, with structural internal loads generally increasing with increasing 
spring stiffness. The final recommended spring stiffnesses are quite high, however, and in 
this range of values the modal and response data are reasonably insensitive to earth 
stiffness modifications. 

The use of beam elements to represent the deck structure is an accurate idealization for 
all deformations except torsion. In the case of deck torsion, because of the torsion- 
bending behavior of the girders, a beam representation is necessarily approximate. The 
nature of ‘■he torsi on -bending action is such that the effective torsional stiffness of the 
deck depends on the torsional mode shape, or wave length, to which the deck is sub- 
jected. This situation makes it possible to determine the deck torsional stiffness with 
acceptable accuracy by calculating the stiffness to correspond to the deck torsional vi- 
bration mode of greatest interest. The deck torsional modes are important principally 
because of their possible involvement in unstable aerodynamic motions (flutter). There- 
fore, the deck torsional stiffness was chosen specifically to obtain accurate modal data 
for the lowest (most flutter-critical) deck torsional mode. The half wave length (one 
lobe) of this mode (see Figure 1) is about 200 feet. Using this length, and postulating 
reasonable girder bending deformations in participation with deck torsional deflections, 
the girder torsion -bending contributions to the deck effective torsional stiffness were 
determined. These contributions are summed with the true torsional stiffness contribu- 
tion, i.e,, those of the slab and the closed box stiffness of the fascia girders, to ob- 
tain the total approximate deck torsional stiffness. 

This procedure necessarily leaves higher deck torsional modes with less accurate (too 
low) torsional stiffnesses, and in general leaves overall bridge modes somewhat in 
error. These errors are negligible since, in the former case, only the lowest deck 
torsion mode was found to be a possible flutter candidate, and in the latter case the 
overall bridge modes are dominated by tower and deck bending and gross deck transla- 
tional influences. It should be noted again that the torsion approximations were neces- 
sitated by the need for a simple, rapidly computed model. To completely resolve the 
deck torsion problem, as was done in the 3D model, would have sacrificed the utility 
of the stick model in the rapid turnaround design support activity. This was an unsatis- 
factory alternative. Moreover, calculations of the stick and 3D model modes confirmed 
the accuracy of the approach used. 

Coding Detoils 

The SAIL (Structural Analyzer Input Language) input language (Reference 2) was used in 
conjunction with NASTRAN bufic data in setting up the structural idealization. Details 
of the coding are described briefly below. 

(i) tbis adjustment was provided by the firm of Shannon and Wilson. 
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Geometry: In using SAIL, the bridge geometry was programmed in the same form 
as it was provided by the design engineers. The grid point coordinates were 
coded in terms of "nose stations", the independent coordinate employed to mea- 
sure distances along the deck cen' .line. The input was greatly simplified through 
the use of NASTRAN's multiple coordinate systems. Cylindrical coordinates were 
used for the portion of the bridge to the left of the tower in Figures 5 and 6, 
which is a circular arc in planview. The remaining straight segment of the bridge 
was input in the rectangular cartesian coordinate system. In the vertical plane, 
the deck describes a parabolic arc of the form 


Z = 155.7 


(Nose Sta - 17070.)^ 
335507 


which was coded directly into the SAIL input deck. An important advantage of 
the SAIL input lies in the fact that variable gridpoint locations and variable num- 
bers of elements are handled in so simple fashion that generation of multiple ideal- 
izations is a minor task. 


2. Multiple Point Constraints: MPC equations were found to be a convenient and 

powerful tool in representing the various connectivities encountered in the bridge 

structure. Structural idealizations using MPC equations are described briefly 

below. 

(a) The cables are rigidly attached to the deck at offset nodes. 

(b) The deck is attached to the tower structure in such a manner that all de- 
grees of freedom except deck vertical bending rotation are required to be 
compatible. In addition, the deck elastic axis is offset (vertically above) 
its supporting cross member in the tower structure, because of the depths 
of the girders, the cross member, and the bearing fitting hardware. 

(c) Similar to the tower attachment described above, the vertically offset deck 
attachment to the piers was enforced by MPC equations. In this case the 
connectivity between the longitudinal motions of the deck and the piers 
was in some designs pinned and in some designs represented by a roller 
support. 

(d) All footings ore connected to the upper ends of the piles by full fixity con- 
ditions enforced by MPC equations. 

(e) In the complex (3D) deck idealization discussed later, MPC equations pro- 
vided the connectivity between the girder webs and the deck slab represen- 
tation. 

The repetitive nature of MPC equations suggests their generation by a subroutine. 

This is discussed briefly under item 3., below. 

3. External Data Generators (EDG): 

(a) The external data generator is a feature within SAIL which provides a 
subroutine type of input generation capability. It is most conveniently 
used for multiple generations of large groups of similar elements and/or 
grid points. In the present problem, this situation occurs for the pile 
foundations. The pile group, including the footing, is a set of 25 elements, 
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29 grid points, plus multipoint and single point constraint conditions. All 
of these input data are prepared by the EDG, in the manner of a subrou- 
tine, needing only onf. set of coding for any number of foundation designs 
to be generated. The use of the EDG permits simple and rapid parameter 
studies on items such as stiffnesses, dimensions, etc., of the pile founda- 
tions. 

(b) The extensive use of multiple point constraints was simplified by creating 
EDG's specifically for the generation of MPC equations. This was done 
for both cylindrical and rectangular cartesian coordinate systems. The 
parameter set for the EDG consists of a list of the two or more nodes to 
be constrained. The EDG recovers the coordinates corresponding to these 
nodes and automatically calculates the constraint equations for a full six 
degree of freedom connectivity. This is particularly useful when cylindri- 
cal equations are employed and in parametric studies where grid point 
changes would otherwise require numerous, potentially erroneous, hand cal- 
culations. The EDG for the pile group generation, described above, calls 
the EDG for MPC equations as required to fix the piles to the footings. 

Computotion Detoils 


The steel bridge idealization shown In Figure 5 was analyzed in 23 different configura- 
tions corresponding to various design changes and parameter studies. The basic model 
consisted of 199 grid points, 95 CBAR elements, and 80 CONROD elements. There 
were 166 MPC equations which in combination with boundary conditions and matrix re- 
ductions reduced the fc'genproblem to 158th order. The runs averaged 6 minutes CPU 
time on the IBM 370 to extract the eigenvalues by Givens' Method and compute 
70 modes. About one-third of the computer time was spent in applying the MPC equa- 
tions. The basic concrete bridge idealization shown in Figure 6 was analyzed in 5 dif- 
ferent configurations. The model consisted of 167 grid points, 84 CBAR elements, and 
64 CONROD elements. There w#'re 149 MPC equations which in combination with 
boundary conditions and matrix reductions reduced the eigenproblem to 14,3rd order. 

The runs averaged 5 minutes 48 seconds CPU time on the IBM 370 to extract the eigen- 
values and compute 70 modes. 
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THREE-D MODEL 



Purpose 

Alfhough the stick model was conceived to be acceptably accurate for both the seismic 
and flutter studies, a refined idealization, the 3D model, was set up for the concrete 
alternate to verify the stick model accuracy. The model was called 3D in reference to 
the idealization of the bridge deck by a slab element (represented as a beam) and indi- 
vidual girder web and flange elements. The beam type idealization of the tower, piers, 
footings, and piles is unchanged from the stick model. A computer plot of the structure 
is shown on Figure 3. 

As discussed earlier, the weakness of the stick model lies in its simulation of deck tor- 
sional stiffness as that of a single member, while in reality the built-up deck resists 
torsion largely through girder bending. Therefore-, the 3D model has as its purpose the 
:accurate representation of girder bending participotion in the overall deck deformations. 

I Description of the Deck Model 

Figure 8 shows schematically three types of behavior of a slab-girder deck. The first 
two apply to a bridge curved in planform, and the last applies for either straight or 
curved decks. All indicate that deck bending, either vertical or horizontal, will couple 
with torsion. The three cases are explained in the text of the figure. Basically, the 
coupling results from two facts: (1) in curved decks, torsion results in lower flange mo- 
tion toward or away from the center of curvature, with a consequent tendency toward 
hoop stresses; (2) in horizontal bending of slab-girder configurations, the elastic shear 
forces are aligned with the shear center of the section (above the deck) while the inertia 
forces are aligned with the mass center. The tendency toward coupling of bending and 
torsion which is described by the figure will affect vibration modes by tending to make 
the mode shapes three-dimensional in character and difficult to identify as pure bending 
or torsional motions. 


In order to represent these coupling tendencies in the finite element mode, it is neces- 
sary to meet several requirements: 

1. Individual girder flanges must be represented in at least axial and horizon- 
tal bending properties. 

2. Diaphragms, cross-bracing, and girder web lateral bending stiffness, all of 
which control lateral motion and therefore hoop forces in the flanges, must 
be modeled. 

3. Flanges must be properly “driven" by the webs; therefore webs must be at- 
tached to the slab in such a way that continuity of displacement and rota- 
tion components is provided. 

4. Structural masses should be properly located. 

All of these requirements were met except the fourth, in the -'location of the masses, 
to simplify the computational problem, the deck moss properties were concentrated at 
the centerline of the slab. For the designs studied, however, the resulting error in mass 
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placement was very small. 

Figure 9 shows the elements used in the 3D model. The deck slab is represented as a 
beam having axial, vertical and horizontal bending and shear, and torsional stiffnesses. 
The centerline of the slab is assigned the deformational freedoms of the deck structure, 
which are the six linear and rotational displacements. All deck motions are constrain 
to these six freeooms. The fascia box girders have a closed cell torsional stiffness. 

This stiffness was added to the torsional stiffness of the deck slab. Flanges are repre- 
sented by beam elements which have axial and horizontal bending and shear stiffnesses. 
The handling of the webs presented difficult problems. It is known, and was further 
verified by calculations, that the NASTRAN plate elements with bending, shear and 
direct stress stiffnesses are of poor accuracy when used as web elements of girders, 
particularly for unsymmetrical cross sections. The erroneous behavior arises from the mei 
brane stiffness of the plate. In order to avoid this difficulty, the girder webs were re- 
presented by combining shear-only plates with bending plates whose only stiffness is 
lateral bending. Because the latter plates cannot maintain the spacing between the 
flanges and the deck, posts are used at the ends of the elements. The axial -force 
stiffnesses (areas) of the girder webs are assigned to the deck slab and to the lower 
girder flanges such that: (1) the elastic axis of the composite deck in vertical bending 
is preserved; (2) the bending moment of inertia of the web of each girder about the 
composite deck elastic axis is preserved. These conditions provide accuracy In girder am 
deck bending and torsional behavior. Axial stretching stiffness of the total deck struc- 
ture, an unimportant factor in the modal analysis, is approximated by these conditions. 

Coding Details 

As with the stick model, the 3D model made use of both external data generators and 
multiple point constraints. The geometry was complicated by the banking of the bridge 
deck (superelevation). Again due to the repetitive nature of the input, SAIL was uni- 
quely suited for data preparation. The principal coding problem is the generation of gr 
point and constraint data for the nine girders. 

The deck centerline geometry and the variable superelevation were computed within the 
SAIL coding, using the equations and data provided by the designers. Using the com- 
puted centerline and superelevation geometrical data in the input parameter set, along 
with component structural data, the deck structure EDG was called. The EDG set up 
the upper girder web (and flange) grid points, the lower girder web (and flange) grid 
points, the girder flange and web elements, the MPC equations which serve to couple 
the girder elements to the six freedoms of the deck centerline, and in addition defined 
the freedoms to be reduced in the eigensolution. The MPC equations rigidly connect 
the upper girder web grid point freedoms to the six freedoms of the grid points on the ' 
deck slab centerline. The EDG is called once for each nose station at which a deck 
grid point is located, thus significantly reducing the magnitude of the coding task. 

This idealization in effect imposes deck cross-sectional bracing (diaphragms) at each dec 
grid point. This Is a correct requirement since the designed diaphragms are located at 
approximately the same nose station spacing as ore the deck grid points. 
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Computotionol Detoils 

The 3D concrete bridge model shown in Figure 3 was subjected to modal analysis. The 
finite element idealization consisted of 623 grid points, 242 CBAR elements, 280 CON- 
ROD element^ 207 CQUADl plates, and 207 CSHEAR webs. There were 1945 MPC 
equations which in combination vi th boundary conditions and matrix reductions reduced 
the eigenproblem to 149th order. A CPU run-time of 19 minutes and 20 seconds on 
the IBM 370 computer was required to extract the eigenvalues by Givens' Method and 
to compute 20 modes. 


SEISMIC AND FLUTTER ANALYSES 

The seismic analysis was performed by the response spectrum method. The full details of 
this method are outside the scope of this paper. Portions of the overall methodology are 
described in Reference 3. The earthquake input data used are in the fom of response 
spectra, and are specifically derived for the West Seattle site conditions . 

The bridge response was determite d in terms of its normal vibration modes. The response 
spectrum method provides maximum individual modal responses to the earthquoke excita- 
tion. Modal summation Is required over very few modes, for most earthquake analyses, 
and is done either as an absolute value sum or a root -square -sum, based on judgement 
and recommendations from past experience (Ref. 3) 

The response spectrum method uses for input the modal analysis data, consisting of vibra- 
tion mode period, generalized mass, and mode shape. In particular, modal response de- 
pends on the degree of coupling between the mode and the uniform vector field which 
describes the motions of the earthquake. This aspect of the seismic analysis requires the 
accessing and processing of very large amounts of structural and modal data. A new 
computer program, used as a NASTRAN post processor, was written to perform the work. 
This program obtains all needed data from the NASTRAN checkpoint/restart tape. The 
set of data reod from the tape consists of files EQEX '1, GPDT, MGG, LAMA, PHIG, 
and OEFl . The complete modal and seismic analysis “xin be done in a single computer 
run, or the NASTRAN and seismic runs can be done separately. The seismic post pro- 
cessing program was found to be very convenient arxl provided a rapid analysis tool. 
Overnight turnaround on combined modal and seismic onolyses was routinely obtained. 

Flutter analysis was done for two types of flutter mechanisms: (1) single degree of free- 
dom stall flutter; and (2) classical bending -torsi on flutter. The calculations were done 
by existing Boeing flutter analysis programs, based on theoretical methods which are be- 
yond the present scope. Aerodynamic data were obtained from wind tunnel tests on 
rroJels of the various bridge deck sections, and modal data were obtained from the NAS- 
TRAN analyses. A subroutine within the seismic program was used to read the normal 
modes from the NASTRAN checkpoint/restart tape and to punch out on cards the required 
rotation and vertical translation displacements of the deck. 

T Data provided by the firm of Shannon and Wilson 
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TYPICAL NUMERICAL RESULTS 


Seismic Analysis 

Each modal analysis computer run provided complete modal data, including element in' 
ternal loads, and also SC4020 plots of the mode shapes. Figures 10 and 11 are com- 
puter plots of the modes which were predicted to be the most important seismic motions 
of the V/est Seattle bridge. The plots shown are for the concrete alternate, but all 
configurations show essentially the same principal types of motions. The first is a lat- 
eral swaying and the second a combined longitudinal - vertical motion which is strongly 
influenced by the cables. The latter is the bridge fundamental mode. These modes 
are important 'leismically for two reasons: (1) their modal frequencies lie in a range of 
strong seismic input; (2) their mode shapes involve essentially unidirectiono! motions of 
the major bridge mqsses, thus detaining strong coupling with the uniform seismic exci- 
tation. 

Figures 12 and 13 show tower moments and shears which were computed for the initial 
steel and concrete designs. The moments and shears shown are those resisting a lateral 
swaying motion, and are caused primarily by modes of the type of Figure 11. The 
stresses for the concrete alternate are larger than those for the steel due mainly to the 
greater deck mass of the concrete design and the close proximity of the concrete modal 
period to a period of strong seismic excitation. These results proved excessively severe 
for strength design purposes. 

As described earlier, the dynamic analysis was continued in support of trade studies for 
the development of a design configuration which is satisfactory for earthquake conditions 
Figure 14 shows the results obtained for a set of seven designs which differ from one 
another primarily in tower configuration. All designs utilize a steel deck structurp. 

From these results tower alternate A was chosen as the recommended configuration . 

The figure lists the modes and modal periods which are critical for both longitudinal 
and transverse earthquake excitations, and gives the resulting maximum tower bending 
moments. The curve shown in the lower right corner of the figure is the earthquake 
response spectrum used in the calculations. 

As discussed earlier- a matter of concern was the effect of the approximation of the 
deck torsional stiffness on the accuracy of the stick model modes, it was for this rea- 
son that the 3D model was used to compute a more accurate set of modal data. Figures 
15 and 16 show 3D model modes corresponding to the stick model modes of Figures 10 
and 11. The agreement in mode shape is excellent. Figure 17 shows the lowest 3D 
model mode in which deck torsion is important. This mode shape justifies the manner 
of computation (the choice of wave length) of the stick model deck torsional stiffness 
which was described earlier. Figure 18 shows a comparison of stick model and 3D 
mcdel modal data for the first ten modes for the concrete alternate. The frequencies 
are tabulated together with a brief description of the modal motions. Note that in 
several coses corresponding modes hove changed order slightly, due to small changes in 
closely spaced frequencies. A careful study of all modal data has shown that in the 

) the defining of these configurations and the choice of tower alternate A were 
done by the firm of Knoerle, Bender, Stone, and Associates, Inc. 
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REPRODUCIBII.ITY OF THE ORIGINAL PAGE IS POOR. 


ef of fhe firsf fen modes only fhe ninfh sfick model mode fails fo agree closely with a 
;orresponding 3D model mode. All other modes show good agreement in both mode 
hope and frequency. 

" •'lutter Analysis 

'he flutter wind speeds weie determined for the initial concrete and steel designs and 
or the current configuration which has been optimized for dynamic response conditions, 
■or the Initial concrete and steel designs, respectively, single degree of freedom stall 
lutter was predicted at steady horizontal windspeeds of 77 miles/hour in the lowest 
orsion mode and 46 miles per hour In the fundamental vertical bending mode. For the 
>ptimized design a torsion stall flutter speed of 244 miles per hour was predicted, with 
he improvement primarily a result of improved aerodynamic shape of the deck section 
ind increased torsional stiffness of the closed box girder design. 

-ONCLUSION 

'he results of the dynamic analyses showed that the initial bridge designs were deficient 
n their ability to withstand a major earthquake or a sustained high wind condition, 
‘hrough the early dynamic analysis parameter studies, however, the directions required 
'or fruitful design modification were defined. A continuing program is in progress to 
mplement these modifications into the design. This woik has resulted in a bridge con- 
Tguration which is satisfactory in resistance to both seismic and wind-induced motions, 
urrently, further design trade studies in conjunction with dynamic analysis are under- 
ay to optimize the design of the lower portion of the tower and the bridge foundation 
or improved earthquake resistance. 

trhis use of the NASTRAN system in the field of civil engineering structures has demon- 
strated a potential for such applications. The benefits of such sophisticated analysis 
pppear particularly great in consideration of the complex structural configurations and 
levere design conditions which are becoming increasingly common in the field of civil 
|sngineering structural design. 
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Flaure la: ELEVATION VIEW 



DECK SUPPORT BEAM 


Figure lb: TOWER CONFIGURATION- ELEVATION VIEW 







Flqure 2; SIMPLE MODEL-CONCRETE ALTERNATE 



Fiqur* 3 : THREE - 0 MODEL — CONCRETE ALTERNATE 
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Fiaure 4a: DECK CROSS-SECTION 




WALL CONSTRUCTION 
BELOW DECK LEVEL 


Fiaure 4b: TOWER CONFIGURATION 


Figure 4: STRUCTURAL CONFIGURATIONS TO RESIST EARTHQUAKE 
AND HIGH WIND CONDITIONS 
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Figure 6 : CONCRETE FINITE ELEMENT (STICK) MODEL 
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Figur* 7 : FINITE ELEMENT PILE ARRANGEMENT 
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Basic Mode of 
Deformation 


Coupling Mode of 
Deformatton 



Lower Flange 
Compression 


lorsioTKil Deformation Causes 
Axial (Hoop) Stress in Flanges, 
Tending to Cause Vertical 
Bending Deformation For 
Stress Relief 



Downward 

Displacement 




Lower Flange 
Tension 


Vertical Displacement 
Causes Axial Stress In 
Flanges, Tending to 
Cause Horizontal and 
Torsional Deformation 
for Stress Relief 



and 

Torsional 

Displacement 


Line of Elastic^ 
Shear Forces 



hertiol Horizontal 
Forces Q^formotlon And 
Forces Involved 


Horizontal Displacement 
Involves Horizontal Shear 
Forces of Elastic Behavior 
and Also Inertial 
(Moss) Behavior - Torsional 
and Horizontal Motions 
Tend to Couple 



Coupled Torsional And 
Horizontal Motion 


Figure 8 : BENDING-TORSION CCXJPLING OF SLAB - GIRDER DECK 
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Deck = Reference Point 

for Displacements 
and Rotations 


\ 


\ 

1 


Element With 
Lateral Bending 
Stiffness Only 


NOTES; 


o Girder Web Axial Area Assigned to Deck and Flange, 
o Box Torsional Stiffneu of Fascia Girder Assigned to Deck. 


0 Deck - Flange Spacing Maintained by Posts at Ends of Elements. 



Figure 9 t ElfMENTS USED IN 3D MODEL 











Fiqure 14: RESULTS OF EARTHQUAKE ANALYSIS TRADE STUDIES 





















































Fiqure 17: LONGITUDINAL, VERTICAL AND TORSION MOTIONS 
- CONCRETE ALTERNATE 


CONCRETE ALTERNATE - CONCRETE ALTERNATE - 

STICK MODEL 3D MODEL 


FREQUENCY (HZ) 

MODE TYPE 

FREQUENCY (HZ) 

MODE TYPE 

,4689 

Y, Z 

.44 

Y. Z 

.7132 

Y. (P) 

.70 

X 

.7305 

X 

.71 

Y. (P) 

.8109 

Y, (P) 

.81 

Y. (P) 

.8537 

Y. Z 

.81 

Y. Z 

.9429 

Y, Z 

.87 

Y. Z 

.9513 

Y. (P) 

.95 

Y, (P) 

1.014 

Y, Z, T 

1.03 

Y. Z, T 

1.065 

X, Z, T 

1.15 

X. Y, Z, T 

1.244 

X, T 

1.19 

X. T 


X - LATERAL 
Y - LONGITUDINAL 
Z - VERTICAL 
T « TORSION 
(P) ■ PIER ONLY 

Figure 18: COMPARISON OF 30 AND STICK MODEL MODAL DATA 
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